Resonant hot charge-transfer excitations in fullerene-porphyrin complexes: a 

many-body Bethe-Salpeter study 
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We study within the many-body Green's function GW and Bethe-Salpeter approaches the neutral 
singlet excitations of the zinc-tetraphenylporphyrin and C7o-fullerene donor-acceptor complex. The 
lowest transition is a charge-transfer excitation between the donor and the acceptor with an energy 
in excellent agreement with recent constrained density functional theory calculations. Beyond the 
lowest charge-transfer state, of which the energy can be determined with simple electrostatic models 
that we validate, the Bethe-Salpeter approach provides the full excitation spectrum. We evidence 
the existence of hot electron-hole states which are resonant in energy with the lowest donor in- 
tramolecular excitation and show an hybrid intramolecular and charge-transfer character, favouring 
the transition towards charge separation. These findings support the recently proposed scenario for 
charge separation at donor-acceptor interfaces through delocalized hot charge-transfer states. 

PACS numbers: 78.40.Me , 78.40.Ri , 71.35.Cc, 73.22.-f 



CO 



-a 
c 

o 
o 



> 
o 
m 

o 



x 



I. INTRODUCTION 

Photo-induced charge-transfer excitations are at the 
heart of a large variety of physical phenomena and re- 
lated applications, from photovoltaic cells for renewable 
energy applications^, to light-harvesting biological sys- 
tems involved e.g. in the photosynthesis processes oc- 
curing in plants or bacterial. In the former case of the 
so-called third generation solar cells such as all-organic 
devices or dye-sensitized Gratzel cells (DSSCs), the large 
binding energy of the photoexcited electron-hole pair re- 
quires the presence of a donor-acceptor interface in order 
to dissociate the bound excitons and produce the free 
carriers at the origin of the photocurrenti. Even though 
the mechanisms explaining such a charge separation are 
still debated^, it is believed that important intermediate 
states are the so-called charge-transfer excitations, with 
the hole and electron located respectively on adjacent 
donor and acceptor units (organic molecules, inorganic 
surface or cluster, etc.) Understanding the structural 
and electronic factors that control the transition from 
the photoinduccd intramolecular (Frenkel) excitons, gen- 
erated by light onto the donor molecules, to a bound 
charge-transfer (CT) state and subsequently to an un- 
bound electron-hole charge-separated (CS) state, is an 
important goal on the road of improving the quantum 
efficiency of solar cells. 

From a theoretical point of view, the correct descrip- 
tion of non-local excitations, such as charge-transfer ex- 
cited states with electrons and holes spatially separated, 
stands as a challenging issue. In particular, it is now 
well established that "standard" time-dependent den- 
sity functional theory (TDDFT) approaches^ based on 
(semi)local exchange-correlation kernels, fail in provid- 
ing a correct description of such charge-transfer (CT) 
excitations.^ This has been demonstrated in the case 
of inter- and intra-molecular CT excitations, but also in 
the case of extended bulk Wannier excitons in semicon- 



ductors for which holes and electrons are separated by a 
large effective bohr radius^ Solutions involving non-local 
kernels based on e.g. range-separated hybrids are being 
actively developed and tested**^ 

In a recent studyii, the lowest charge-transfer exci- 
tation in the zinc tetraphcnylporphyrin (ZnTPP) donor 
associated with the C70 fullerene acceptor was studied us- 
ing a density functional theory (DFT) approach within 
an elegant and computationally efficient constrained for- 
malism with a penalty functional forcing the jump of a 
chosen fraction of an electron from the ZnTPP highest 
occupied molecular orbital (HOMO) to the C70 lowest 
unoccupied molecular orbital (LUMO). It was shown in 
particular that the related excitation energy of ~ 1.8 eV 
in case of a full electron transfer was much larger than the 
~ 0.9-1.0 eV energy obtained with TDDFT based on the 
PBE kernel^ Even though the available experiment^ 
was performed in a solvant with a different porphyrin 
derivative, the constraincd-DFT result is in much better 
agreement with the range of experimental values. 

Computationally very efficient and accurate, such 
DFT-based techniques^ cannot provide however the de- 
tailed excitonic (or absorption) spectrum at the donor- 
acceptor interface. In particular, the energy position 
of the donor intramolecular (Frenkel) low-lying exciton 
with respect to the charge-transfer excitonic spectrum, 
including the higher energy (hot) charge-transfer states, 
is a crucial information needed to understand the dis- 
sociation of the excitons generated on the donor side 
and that have diffused to the donor-acceptor interface. 
Mechanisms where intramolecular excitons evolve into 
charge-separated states through hot (high energy) CT 
states have been recently demonstrated experimentally^ 
so that the knowledge of the entire excitation spectrum 
(at least in a few electronvolts energy window) may be of 
crucial importance. Besides the energy spectrum alone, 
the related knowledge of the excitonic wavefunctions is 
a crucial ingredient governing e.g. the calculation of 
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charge-transfer rates. Furthermore, in the case of the 
constrained DFT study of the C7o-ZnTPP complex, the 
underlying assumption of a full electron transfer may de- 
serve closer inspection. 

Recently,— ~— charge-transfer excitations in small gaz 
phase donor-acceptor complexes, associating tetracya- 
noethylene with several acene derivatives, were stud- 
ied using the many-body perturbation theory frame- 
work within the so-called GW— ~— and Bcthc-Salpctcr 
(BSE) 2 ^— formalisms. It was shown in particular— ii£ 
that an excellent agreement with experiment 2 ^ could 
be obtained for the lowest charge-transfer excitation, 
with a mean absolute error within 0.1 eV. This accu- 
racy matches that of the latest range-separated func- 
tional TDDFT studies with the range parameter ad- 
justed to provide an accurate quasiparticle band gapi^ - 
The same formalism was used to study intramolecular 
charge-transfer excitations in a recent family of molec- 
ular dies for DSSCs, the coumarins 3 - 1 -, showing a mean 
absolute error of the order of 0.06 eV as compared to 
reference quantum chemistry coupled-cluster— or tuned 
range-separated hybrid TDDFT calculations^ 3 ^!. 

In the present work, we report the detailed excitonic 
spectrum of the C70 fullerene associated with the zinc 
tctraphcnylporphyrin (ZnTPP). This study is performed 
within the many-body Green's function GW and Bethc 
Salpetcr formalisms using the same cofacial geometry as 
that used in the constrained DFT study of Ref. [Tl|. Our 
results are in excellent agreement with constrained DFT 
calculations for the lowest lying CTo charge-transfer ex- 
citation, but provide further the entire excitation spec- 
trum. It is shown in particular that several hot electron- 
hole excited states lie in between the lowest intramolec- 
ular ZnTPP* donor excitation and the CTo exciplcx. Of 
peculiar importance, several hot states with hybrid in- 
tramolecular and charge-transfer character are found to 
be resonant in energy with the photo-induced ZnTPP* 
exciton, a finding important in the perspective of the re- 
cently proposed scenarios where charge-separation occurs 
directly through delocalized hot charge-transfer states. 



II. METHODOLOGY 

Our calculations are performed with the Fiesta 
codeiS^"— , a recently developed GW-BSE package 
based on a Gaussian auxiliary basis expansion of 
the needed two-body operators such as the dielec- 
tric susceptibility, and the bare or screened Coulomb 
potentials 3 ^. Dynamical correlations are calculated ex- 
actly using contour deformation techniques without any 
plasmon pole approximation. The needed Kohn-Sham 
eigenstates are obtained from the Siesta package 3 ^ 
with a well- converged triple-zeta plus double polar- 
ization basis (TZDP) and standard norm-conserving 
pscudopotcntials^ - combined with the local density ap- 
proximation (LDA)^i for the exchange-correlation func- 
tional. All empty states are included in the summa- 



tion over transitions involved in the calculation of the 
independent-electron susceptibilities. As shown recently 
in the study of charge-transfer excitations in several w- 
conjugated organic donor-acceptor systems^ the present 
scheme provides results in excellent agreement with 
planewave-based GW-BSE studies performed with the 
Yambo codo 4 ^ exploiting Kohn-Sham eigenstates gener- 
ated with the Quantum-Espresso package^ 3 -. The defini- 
tion of the auxiliary basis, consisting of six Gaussians per 
s,p,c?-channel, together with extensive convergency tests 
as a function of the auxiliary basis size, can be found 
in a recent study of a photosynthetic donor-acceptor 
complex^. 

We go beyond the standard "single-shot" Go Wo calcu- 
lation by performing a partial self-consistent loop onto 
the quasiparticle eigenvalues, leaving the Kohn-Sham 
wavefunctions unchanged. Such a scheme has been 
shown for molecular systems to lead to quasiparticle 
energy gaps, between the highest occupied (HOMO) 
and lowest unoccupied (LUMO) molecular orbitals, in 
better agreement with experiment than the standard 
GqWq calculation based on DFT-LDA or DFT-PBE 
eigenstate a 35 i 44 ~ — at a cost much reduced as compared 
to full self-consistency4£ Similarly, GW-BSE CT excita- 
tion energies in reference donor-acceptor complexes were 
shown to agree much better with experiment if such a 
partial self-consistency is performed at the initial GW- 
level . 16 ' 17 The dependence of the final GW quasiparticle 
energies on the starting Kohn-Sham eigenstates is sig- 
nificantly reduced, and GW calculations starting from 
DFT-LDA or Hartree-Fock-like 4 ^ eigenvalues were shown 
to converge towards the same GW spectrum for several 
organic molecules . 35 ' 50 

While the GW calculations provide accurate quasipar- 
ticle energies, including not only the correct HOMO- 
LUMO gap, but also a correct spacing and ordering 
within the occupied (unoccupied) states, the subsequent 
resolution of the Bcthe-Salpeter equations (BSE) aims 
at providing the correct neutral excitation energies, ac- 
counting in particular for electron-hole interaction. We 
go beyond the Tamm-Dancoff approximation, including 
namely the coupling of resonant and non-resonant exci- 
tations. As shown in previous work , 31 ' 51 ' 52 the Tamm- 
Dancoff approximation leads in the case of molecular (or 
nanosized) systems to an overcstimation of the transition 
energies by as much as 0.3 eV. 

The BSE Hamiltonian is most commonly expressed in 
the \4>t<i>^ > product basis, where (</>f ) and (<p^) are unoc- 
cupied (electron) and occupied (hole) Kohn-Sham states, 
respectively In particular, the resonant contribution is 
composed of a "non- interacting" diagonal part (H dla9 ), 
a direct (7J dlrect ) anc j an exchange contribution (H exch ) 
term, with: 

llIX = - J drdr>f(r)^(r')W(r,r')^(r)^(r'), 



H 



= 2 / rfrdr'^(r)^(r) V (r,r')0 e fe (r')^(r'), 



where the (eSJ) are the quasiparticle (GW) energies, 
while W(r,r') and v(r,r') are the (statically) screened 
and bare Coulomb potential, respectively. It can be no- 
ticed in particular that the direct term does not cancel 
in the limit of spatially non-overlaping electron and hole 
states, an important remark concerning charge-transfer 
states. 

Due to the fast increase in size with the number of 
unoccupied states of the electron-hole two-body product 
\4>f cj>j > basis states, we restrict our BSE basis to tran- 
sitions between the occupied states and 360 unoccupied 
levels, that is namely empty states up to 20 eV above 
the HOMO level. Convergency tests can be found again 
in Ref. [38| for donor-acceptor systems of similar size. We 
emphasize that the low-lying charge-transfer excitations, 
with much weight on the lowest occupied and unoccu- 
pied Kohn-Sham electron-hole product states, converge 
extremely quickly with the size of the BSE basis, but this 
is not necessarily the case for intramolecular excitations 
or higher (hot) CT excitations. 
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III. RESULTS 

We adopt the same geometry as that used in the 
constrained-DFT study of the C7o-ZnTPP complex in 
Ref. [111. Namely, the distance between the Zn atom 
and the center of the closest electron-rich (6:6) C70 car- 
bon bond is fixed to 2.8 A using a constrained relax- 
ation scheme. Such a geometry stems in particular 
from the analysis of the stability of several fullerene- 
porphyrin complexes in various geometries^. How- 
ever, as shown recently in the case of a fullerenes- 
polythiophcnes interface^ significant disorder is to be 
expected at the bulk donor-acceptor interface, so that 
the present geometry is adopted for sake of comparison 
with available previous calculations of charge-transfer ex- 
citations. 



A. Quasiparticle spectrum within GW 

We show in Fig. Q] a cartography of the Kohn-Sham 
(left) and GW (right) eigenvalues together with a rep- 
resentation of a few associated eigenstates. As ex- 
pected, the quasiparticle HOMO-LUMO gap is signif- 
icantly opened from the DFT-LDA to the GW calcu- 
lation. The quasiparticle GW HOMO-LUMO gap is 
found to be of 3.45 eV, namely 2.6 eV larger than the 
corresponding DFT-LDA Kohn-Sham result. Our GW 
HOMO-LUMO gap is in very good agreement with the 
3.56 eV value found for the same system in Ref. [H by 
taking differences of total energy between the neutral and 
positively or negatively charged complex at the DFT- 
PBE level within the standard ASCF approach. Simi- 



FIG. 1: (Color online) Cartography of the Kohn-Sham (left) 
and GW (right) eigenvalues. The energy levels are distributed 
into a (donor) or (acceptor) column according to the spatial 
localization of the corresponding eigenstates. Frontier Kohn- 
Sham orbitals are represented next to their associated energy 
level. 



larly, our ionization potential of 6.25 eV is in close agree- 
ment with the reported 6.21 eV ASCF-PBE value. In 
the absence of gas phase experimental data, such a close 
agreement between the two approaches is very satisfac- 
tory. 

The donor-acceptor, or type-II, character of the com- 
plex is evidenced by the relative positions (frontier or- 
bitals energy offsets) of the C70 and ZnTPP eigen- 
states across the gap. At the GW level, the ZnTPP 
HOMO (LUMO) is located 0.95 eV (1.05 eV) above 
the C 70 HOMO (LUMO) stated The (HOMO) and 
(HOMO-1) orbitals are the standard porphyrin Gouter- 
man orbitals^, the (HOMO-2) being also localized onto 
ZnTPP, while the (nearly) 3-fold level beneath derive 
from the HOMO of the C 70 fullerene^^ The three low- 
est unoccupied levels, (LUMO) to (LUMO+2), derive 
from the 2-fold e\ and a\ unoccupied levels of C70, the 
two unoccupied Gouterman orbitals located on ZnTPP 
being the (LUMO+6) and (LUMO+7) levels. Such a 
mapping of orbitals will allow us to analyze the optical 
excitations in terms of intramolecular (Frenkel) versus 
charge-transfer transitions. 

The intramolecular ZnTPP HOMO-(LUMO+6) gap 
is found to be 4.5 eV, of comparable size with the in- 
tramolecular C 70 (HOMO-3)-LUMO gap (4.4 eV), but 
~1.1 eV larger than the true charge-transfer HOMO- 
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LUMO quasiparticlc gap. The relative strength of the in- 
tramolecular (Frenkel) and charge-transfer electron-hole 
pairs binding energy will reshuffle the ordering of the 
photo-induced transition as compared to the GW occu- 
pied to unoccupied quasiparticle energy differences. The 
comparison between the ~ 3.5 eV HOMO-LUMO gap 
and the 1.6-1 .8 cV optical absorption threshold found in 
Refs. [TTIIl3l indicates that in such systems the electron- 
hole interaction is extremely large, amounting to about 
half the quasiparticle gap. This is what we now study 
within the Bethe-Salpeter formalism that aims at ac- 
curately describing the electron-hole interaction in the 
photo-induced excited states. 



B. The low-lying charge-transfer state within BSE 

Our GIT-BSE results for the neutral (singlet) excita- 
tion energies are provided in Fig. [2]where we have distin- 
guished the excitations with respect to their intramolecu- 
lar or charge-transfer character. Such a distinction stems 
from the analysis of the weight of the Bcthc-Salpctcr 
two-body eigenstates onto the Kohn-Sham occupied- 
unoccupied product state basis. We provide further, in 
Fig. [2] and Fig. [3l a representation of selected two-body 
excitonic eigenstates ipf, h (r e ,rh), with (r e ,Th) the elec- 
tron and hole positions, respectively, and v the excita- 
tion index, by plotting in grey the hole-averaged electron- 
distribution, that can be obtained by taking the expecta- 
tion value of the electron S(r — r e ) position operator over 
tpf, ■ Similarly, the electron- averaged hole distribution is 
represented in blue (wireframe). 

The G1T-BSE low-lying singlet excitation is a clear 
charge-transfer state located at 1.65 eV above the 
ground-state (we will use the notation CT$ in the follow- 
ing). Since the C70 LUMO is (nearly) 3- fold degenerate 
(see Fig. 1), the 3 low lying CT excitations are linear 
combinations of transitions between the HOMO and the 
LUMO, LUMO+1 and LUMO+2 levels, namely: 



1,2,3 



(r e ,r fe ) 



n=0 



(re). 



For the lowest lying transition, we find typically: (\c\\ 2 = 
0.27, \c\\ 2 = 0.42). As such, these low-lying 



,1|2 



0.29, 

excitations arc very clear charge-transfer states with 
the jump of nearly one electron from the donor related 
HOMO level to the manifold of the acceptor LUMO. This 
reflects in the representation of the two-body wavefunc- 
tion for the lowest CTq state (Fig. UK) showing that the 
(hole-averaged) electron density is very clearly located 
onto the C70 molecule while the (electron-averaged) hole 
density resides onto the donor ZnTPP molecule. 

Our ^1.65 eV G1P-BSE excitation energy for the low- 
est excitation is in good agreement with the 1.8 eV value 
extrapolated from constrained-DFT calculation^ under 
the assumption that a full electron is transfcred between 
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FIG. 2: (Color online) Cartography of the BSE singlet neu- 
tral excitation energies. The transitions are symbolically lo- 
cated from the left to the right in terms of their donor in- 
tramolecular, charge transfer (CT) or acceptor intramolecular 
character. In particular, transitions with a representative en- 
ergy level displaced from the center to the left (right) indicates 
charge-transfer excitations with some weight onto intramolec- 
ular ZnTPP (C70) transitions (see text). A few selected two- 
body eigenstates are represented (see text and Fig. [3}caption 
for details). 



the donor and the acceptor. On general grounds, apart 
from the performances of the PBE functional in cal- 
culating differences of total energies for such systems, 
a constrained approach within a variational framework 
necessarily yields an upperbound for the excitation en- 
ergy. As emphasized further in Ref. [Til the assumption 
of a slightly reduced electron transfer would bring the 
constrained-DFT and BSE calculations in even better 
agreement. An analysis of Fig. 5 of Ref. [ll| shows that 
the largest charge transfer explicitely studied by Ghosh 
and Gebauer corresponds to a transfer of ^0.97 electron 
with a corresponding excitation energy of ~1.7 cV, in 
close agreement with our findings, the 1.8 eV value being 
an extrapolated energy for full electron transfer. 



As observed in Refs. lllil55l the direct comparison with 
experiment is difficult. Existing experimental values by 
Mukherjee and coworkers^ report the lowest excitation 
for the parent octadecyloxyphenyl-H2TPP (freebase por- 
phyrin) G70 at 1.7 eV in toluene solvant and 1.8 eV in 
chlorophorme. As shown in Refs. loSllaOl the replacement 
of the freebase porphyrin by a central Zn atom leads 
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FIG. 3: (Color online) Symbolic representation of the 
ip^, h {v e ,rh) two-body excitonic BSE eigenstates (see text), 
with (a) v = (CTo), (b) v = 15 (lowest intramolecular 
ZnTPP* excitation), and (c,d) v = 16 and v = 17 hot CT n 
excitations with both a charge-transfer and an intramolecular 
character. The grey isocontours represent the hole-averaged 
electron distribution of charge and the wireframed blue con- 
tours the electron-averaged hole distribution. 



to a 0.1-0.2 eV redshift of the CTo excitation energy. 
On the contrary, removing the effect of solvant would 
lead to a compensating 0.2-0.3 eV blueshift (see analysis 
in Ref. [poh . even though the redshift of the CT excita- 
tion energy observed by Mukherjee and coworkers upon 
changing chlorophorme by the less polar toluene solvant 
indicates potentially more complex trends. 

Even though the comparison with experiment is diffi- 
cult, we oberve that the agreement is significantly better 
than what is obtained within TDDFT and (semi)local 
kernels, namely ~0.9 eV in the TD-PBE study of Ref.GH. 
The underestimation of charge-transfer excitations en- 
ergy within TDDFT and (semi)local kernels is now well 
documented. Such a discrepancy lies in the lack of 
electron-hole interaction in the case of non-overlaping 
HOMO-LUMO states, a situation that can be clearly evi- 
denced e.g. by the absence of dispersion of the CT excita- 
tion energy as a function of donor-acceptor distance^^ 
contradicting the simple electrostatic limit of a (-1/D) 
scaling of the interaction between two well separated 
positive and negative charges, where D is some aver- 
aged measure of the electron-hole distance. The lack of 
electron-hole interaction for CT states indicates that the 
CTo excitation energy reduces to the (too small) Kohn- 
Sham HOMO-LUMO gap, that is indeed the -0.85 eV 
HOMO-LUMO gap that we find at the DFT-LDA level. 



C. Testing models for charge-transfer excitations 

The very clear charge-transfer character of the CTq 
low- lying excitation allows further to test approximations 
of use in the quantum chemistry community by extract- 
ing from the BSE Hamiltonien the following matrix ele- 
ments of the bare and (statically) screened Coulomb op- 
erator, namely: 



E 



B,(l) 
eh 



E 



B,(2) _ 



drdr'|^ omo (r)| 2 M/(r, r')|tf» lumo-\-n \I J | ? 
drdr' |</>/ lomo (r)| 2 i;(r,r')|(/) ; 

um o-\-n \* J | 5 



with (n=0,l,2). Such matrix elements are nearly in- 
dependent of (n) with: £ e B h ' (1) =1.56 ± 0.02 eV and 

£^' (2) =1.72 eV ± 0.06 eV. The exchange matrix ele- 
ments, namely: 



E eh 



drdr'(j) h (r)<j)i(r)v(r,r')<t> h (r')(l>i(r') 



where (h) stands for (homo) and (1) for either (lumo), 
(lumo+1) or (lumo+2), amount only to less than 1 meV 
as a result of the charge separation between the (HOMO) 
and the (LUMO) states. Subtracting now the screened 

B (1) 

Coulomb electron- hole interaction E el l to the GW 
quasiparticlc HOMO-LUMO gap (3.45 eV), one obtains 
a simple estimation of the CT excitation energy, namely 
1.89 eV, in reasonnable agreement with the 1.65 eV value 
obtained with the full GVF-BSE treatment. As an even 
simpler approximation, which does not necessitate the 
calculation of the full dielectric matrix, subtracting now 

B (2) 

the bare Coulomb potential related E e ^ energy to the 
GW HOMO-LUMO gap, one obtains 1.72 eV in (acciden- 
tally) even better agreement with the GTU-BSE value. 

The present results clearly suggest that even in the 
limit of a separation of —2.8 A in a cofacial geometry, 
namely a distance and a configuration maximizing the 
wavcfunction overlap between the donor and the accep- 
tor for such 7r-stacking systems, the rather crude "elec- 
trostatic" model consisting in replacing the electron-hole 

B (2) 

binding energy by the bare Coulomb integral E ,' v pro- 
vides very acceptable results, with an error as compared 
to the full BSE value of the same order of magnitude than 
the difference between the GTU-BSE treatment and the 
(extrapolated) constrained-DFT approach. This simple 
evaluation of the electron-hole binding energy in charge- 
transfer states was used for example in Ref. [sTj to pre- 
dict the evolution of the electron-hole binding energy for 
several donor-acceptor complexes and relative geometry. 
This simple scheme relies however on the assumption that 
the excited state is a clear charge-transfer state, with no 
mixing with higher energy (0f<^) single-particle transi- 
tions, an approximation of the same nature than the one 
used in the constrained-DFT approach with the "guess" 
that a full electron is transferred. 
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D. "Hot" excited states: resonant intramolecular 
and charge-transfer excitations 

While our GW^-BSE calculations yield results in 
good agreement with the computationally more efficient 
constraincd-DFT approach for the low-lying CTq level, 
the GW^-BSE framework provides the entire excitation 
spectrum as illustrated in Fig. 2. In particular, the 
relative position of the lowest intramolecular (Frenkel) 
donor excitation (labeled ZnTPP* here below), and the 
charge-transfer excitations at the donor-acceptor inter- 
face, is believed to be a crucial factor in understand- 
ing the mechanism(s) controlling the transition from a 
strongly bound electron-hole pair to a charge-separated 
(CS) state with unbound holes and electrons diffusing to 
the electrodes. In recent combined experimental and the- 
oretical studies^i£2, the role of hot CT„ charge-transfer 
excitations in mediating the transition between Frenkel 
and CS states was advocated on the basis of experimental 
pump-probe experiments and quantum chemistry calcu- 
lations. Located at higher an energy than the lowest ly- 
ing CTq level studied here above, the CT n hot states are 
significantly more delocalized due to their extra kinetic 
energy, a favourable feature for evolving into a charge- 
separated state£i£2 

We can now analyze the present ZnTPP-C7o system 
along these lines and study in particular the position 
of the ZnTPP* lowest intramolecular excitation with re- 
spect to the manifold of CT states. Within the present 
GW-BSE approach, the ZnTPP* Q-like^ excited state is 
located at ~ 2.14 eV (see left "donor" column in Fig. ^ 
with a corresponding two-body eigenstate provided in 
Fig. 02h). Such a state is composed of transitions from 
the (HOMO-n) (n=0,l,2) ZnTPP-related levels to the 
(LUMO+6) and (LUMO+7) levels deriving from the low- 
est unoccupied levels of the isolated ZnTPP molecule (see 
Fig. [I]). The GW-BSE excitation energy (~ 2.14 eV) is 
very close to the experimental Q xs value of 2.09 eV for 
gas phase ZnTPPi^ Previous GW-BSE calculations of 
Q-transitions in a parent zincbacteriochlorin molecule^ 
and freebase H2TPP— were shown to lead as well to an 
agreement within 0.1 eV as compared to experiment for 
such intramolecular low-lying excitations. 

As depicted in Fig. [5J the ZnTPP* intramolecular ex- 
citation lies near or above several charge-transfer excita- 
tions located above the CTo state. Among these higher- 
energy, or hot, excited states, interesting transitions are 
the excitations (y = 16) and [y = 17) that lies within 
10 meV of the ZnTPP* lowest intramolecular excitation. 
With the notations defined above, one finds for the cor- 
responding two-body eigenstates: 

VV=16,17 ( r e,r/j) ~ c£,J+6 <t>homo{*h)<t>lumo+&{*e) 
^h— 1,2+7 ^homo—1 
~r 2 I Wfiomo—z (r h )^iu m o{r e ), 

with (K ;+6 | 2 ,K_^ +7 | 2 ,K_ 2 /) equal to 



(0.40,0.32,0.15) for (v = 16) and (0.10,0.09,0.74) 
for {v = 17). Again, the (lumo+6) and (lumo+7) 
levels are mainly the nearly degenerate lowest unoc- 
cupied Gouterman orbitals of the ZnTPP molecule. 
As such, the two-body wavefunctions present both a 
charge-transfer character, with the (homo-2) to (lumo) 
transition component, but also an intramolecular ZnTPP 
character, with the (homo) to (lumo+6) and (homo-1) 
to (lumo+7) components. The corresponding excitonic 
wavefunctions are represented in Fig. 3(c,d), evidencing 
this hybrid character. The (y = 16) excitation is mainly 
a ZnTPP intramolecular transition showing however a 
significant probability for the electron (grey contours) to 
jump onto the C70 acceptor. The (y = 17) excitation is 
mainly a CT transition with the excited electron sitting 
mainly onto the C70 acceptor, but with a strong weight 
as well onto the ZnTPP donor, in great contrast with 
the low lying CTo state which is a pure charge-transfer 
state. 

While the present system, comprising only one donor- 
acceptor complex, cannot be used to test the dereal- 
ization of the hot excited states over several donor (for 
the hole) or acceptor (for the electron) molecules, a de- 
localization expected to favour charge separations^ the 
present finding indicates that the extra kinetic energy 
of the hot excited states, as compared to the low ly- 
ing CT state, induces also the delocalisation of the 
hole and electron over the neighbouring donor and ac- 
ceptor molecules, generating hybrid intramolecular and 
charge-transfer states. The present results are consis- 
tent with a previous semi-empirical stud}*^ of a Cqq- 
polymcr blend evidencing a rapid drop of the "ionicity" 
of charge-transfer excited states with increasing energy, 
a reduction of the charge-transfer character characteriz- 
ing an hybridation with intramolecular transitions. It 
was shown however that the number, energy and char- 
acter of the CT states lying below the intramolecular 
donor (polymer) excitations were strongly dependent on 
the parametrization of the intersite (Vij) Coulomb po- 
tential, making the case for accurate parameter-free ab 
initio calculations. 

Such an energy alignment and wavefunction overlap 
between the ZnTPP* intramolecular low- lying excita- 
tion and these hot CT„ states is a strong indication 
that transitions between these resonating states will be 
strongly favoured^ The present findings may suggest 
an additional criteria for designing donor-acceptor com- 
plexes with efficient charge separation through tran- 
sient hot (delocalized and weakly bound) charge-transfer 
states, namely the existence, at the energy of the low- 
est ZnTPP* donor intramolecular excitation, of a large 
density of resonating states with hybrid intramolecular 
and charge-transfer character. Adopting a simple hy- 
drogenoid model for the electron-hole pair, this can be 
expected to be achieved by maximizing the energy dif- 
ference between the low-lying CTq charge-transfer state 
and the intramolecular ZnTPP* donor excitation. Even 
though equivalent to the standard criteria for energy- 



gradient-driven transitions to the CTo state, we empha- 
size that the argument here applies to a scenario where 
electron-hole dissociation occurs through hot charge- 
transfer states, avoiding as much as possible the ad- 
verse recombination in the low lying strongly bound CTo 
exciplex.— ' 62 i 67 Following Ref. 0, we represent schemati- 
cally in Fig.|3]a flow of transitions from donor intramolec- 
ular to resonant hot CT„ states, and finaly to a charge- 
separated (CS) state. 

In a recent GIU-BSE study of a Ceo-quaterthiophene 
system^ much emphasis was given to calculating the 
difference of energy (AO) between the lowest donor 
intramolecular Frenkel (FE) exciton and the charge- 
transfer CT state. AO was found to be sensitive to 
the relative donor-acceptor orientation, with an average 
positive value stronly detrimental to the FE — > CTo 
transition. Further, as observed here above, the bind- 
ing energy (Eg ) of the electron- hole pair in the CTq 
state was found to be much larger than room tempera- 
ture. Even though based on very accurate calculations, 
such results do not allow to understand the observed 
efficiency of state-of-the-art fullerene-polythiophene sys- 
tems. While the difference between isolated gas phase 
donor-acceptor dimers, on which calculations are usually 
performed, and more realistic bulk donor-acceptor het- 
crojunctions (BHJ) was invoked^ a possible alternative 
interpretation lies in the scenario where charge separation 
occurs directly through hot CT„ states, without transit- 
ing to the highly bound CTo state, so that variations in 
AO and Eq TQ become less relevant. 

Back to the present ZnTPP-Cyo system, the further 
analysis of the CT states lying between the CTo and 
ZnTPP* levels shows several states with hybrid charge- 
transfer and intramolecular character, but with a clear 
weight onto intra-fullerene excitations. Since photoex- 
cited pairs are mainly generated in the donor (ZnTPP) 
part of the heterojunction, such excitations cannot be 
efficient intermediate states for charge separation, indi- 
cating on the contrary a pathway to an adverse relocal- 
ization of the electron-hole pair onto the fullerene side. 

As a final analysis, we now compare our results 
to a recent study of the excited states of fullerene- 
porphyrin complexes^ based on another computation- 
ally efficient DFT-based framework. As described in 
Ref. [H, this efficient, but approximate, approach re- 
lies on a second-order perturbative correction to the 
"excited" Slater determinant obtained by replacing one 
occupied single-particle DFT Kohn-Sham orbital by 
an empty one in the ground-state Kohn-Sham many- 
body wavefunction. Since the empty and occupied or- 
bitals considered were the (HOMO) to (HOMO-2) and 
(LUMO) to (LUMO+2) orbitals, only low-lying charge- 
transfer states could be studied, excluding in particu- 
lar the intramolecular ZnTPP excitations, involving the 
(LUMO+6) and (LUMO+7) single-particle empty or- 
bitals, and the hybrid states above-mentioned. 

The lowest charge-transfer CTq excitation within this 
scheme was found to be located at 1.95 eV above the 




FIG. 4: (Color online) Symbolic representation (following 
Ref. 0) showing the potential flow of transitions from in- 
tramolecular (Frenkel) excitations, localized on the donor, to 
resonant hot charge-transfer (CT n ) states (red arrow) and 
subsequently unbound charge-separated (CS) electron-hole 
pair (orange arrow). The adverse decay into the low lying 
CTo state is also represented with the green arrow. 



ground-state.— Our GVF-BSE value is ~ 0.3 eV smaller. 
Even though slightly less satisfactory than the compari- 
son with the variational constrained DFT calculation^, 
the agreement remains reasonable. Since this approach is 
perturbative, it is difficult to conclude whether, on gen- 
eral grounds, it should produce an upper- or lower-bound 
for the exact result. In however excellent agreement with 
this perturbative technique, the three lowest excitations 
are CT excitations with a strong weight on the HOMO 
to LUMO, LUMO+I and LUMO+2 quasidegenerate lev- 
els, the following excitation being mainly a (HOMO-I) to 
(LUMO+2) transition. Inclusion of intramolecular tran- 
sitions in the pool of studied excitations may certainly 
help in favouring the comparison with other techniques 
and in revealing the excitations with hybrid intramolec- 
ular and charge-transfer character. 



IV. CONCLUSION 

We have studied within the many-body Green's func- 
tion GW and Bethe-Salpeter formalisms the singlet ex- 
citation energies in a fullerene-porphyrin C7o-ZnTPP 
complex. The lowest-lying excited state is a charge- 
transfer (CTo) state with a ~ 1.65 eV energy, a value 
in good agreement with recent constrained DFT calcula- 
tions leading, for the same geometry, to a 1.8 eV upper- 
bound in the asymptotic limit of a full electron transfer. 
The energy of the CTo state can be estimated simply 
from the knowledge of the correct GW HOMO-LUMO 
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gap and of simple Coulomb integrals between the HOMO 
and LUMO Kolm-Sham eigenstates density. Beyond the 
lowest excitation, the GW-BSE framework provides the 
full excitation spectrum, allowing in particular to assess 
the relative position of the lowest intramolecular ZnTPP* 
donor excitation with respect to the manifold of hot 
charge-transfer CT n states. We reveal in particular the 
existence of several excitons with an hybrid intramolec- 
ular and charge-transfer character which are resonant in 
energy with the ZnTPP* excitation. Such findings sug- 
gest a rapid transition from the photo-generated ZnTPP* 
state to hot CT n states with excess kinetic energy as com- 
pared to the lowest CTo state, favouring thus the tran- 
sition towards unbound charge-separated states. More 
direct comparison to experimental results requires to ac- 
count for environmental effects, namely the renormalisa- 
tion of the quasiparticle and neutral excitation energies 



from the gas phase to realistic bulk phases. The combi- 
nation of the present GW and Bethe-Salpeter formalisms 
with proper and accurate polarizable models is certainly 
an important goal for the community in order to prop- 
erly discriminate between the various scenario proposed 
for charge separation in a realistic environment. 
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